Accuracy and range of validity of the Wigner surmise for mixed symmetry classes 

in random matrix theory 
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Schierenberg et al. [Phys. Rev. E 85, 061130 (2012)] recently applied the Wigner surmise, i.e., 
substitution of oo x oo matrices by their 2x2 counterparts for the computation of level spacing 
distributions, to random matrix ensembles in transition between two universality classes. I exam- 
ine the accuracy and the range of validity of the surmise for the crossover between the Gaussian 
orthogonal and unitary ensembles by contrasting them with the large-A^ results that I evaluated 
using the Nystrom-type method for the Fredholm determinant. The surmised expression at the 
best-fitting parameter provides a good approximation for < s < 2, i.e., the validity range of the 
original surmise. 
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Recently, Schierenberg et al. [l[ publicized an exhaus- 
tive study on the application of the Wigner surmise [2[ to 
the computation of level spacing distributions (LSDs) of 
random matrix ensembles in "mixed symmetry classes," 
i.e., 

H = Hp + XHfji £ {2 x 2 Hermitian matrices} (1) 

with Hp and Hp* taken from the Gaussian measures 
dHp q~^ 1 / 2 ^ h p of Dyson indices j3 and /?', respectively, 
succeeding the works by Lenz and Haake and others 
Although the substitution of random matrices of infinite 
rank with those of the smallest possible rank appears too 
bold an ansatz, the resultant distributions P{s) of un- 
folded level spacings s = (ei+i — e i)/A [e, G Spec(Hp), 
A: mean level spacing around ej] for the three Dyson 
universality classes (3 = 1 , 2 and 4 arc well known to be 
in good agreement with the exact large- TV results (HQ, 
especially in the range < s < 2 (Fig. 1). 

The intuitive reason for this is that the small- s behav- 
iors of LSDs are dominated by the repulsion of two (oth- 
erwise colliding) eigenvalues of a Hamiltonian under per- 
turbation. Eigenvalues adjacent to the two anticrossing 
eigenvalues are expected to have very little influence on 
this universal repulsive behavior, leading to the validity 









— GOE " 




— GUE 




GSE 







of the 2x2 approximation for relatively small spacing s. 
However, even this intuitive account can hardly be taken 
for granted in the case of Hamiltonians in universality 
crossover. Nevertheless, the authors of [l[ have reported 
that the numerically generated spectrum of random ma- 
trices of large rank N(~ 400), 

H = Hp + aH (i i £ {N x N Hermitian matrices} , (2) 

can be fitted well to the Wigner-surmised expression from 
([1]) whose A parameter is close to A = (sp/A)a (Figs. 5, 
7, 9, and 10 of Here, sp is the mean level spacing of 
the unperturbed 2x2 random matrix Hp. 

On the other hand, I have recently evaluated LSDs for 
crossovers between Gaussian orthogonal (GOE) and uni- 
tary (GUE) ensembles and between Gaussian symplectic 
ensembles (GSE) and GUEs {(3 = 1,4, (3 1 = 2) in the 
large- N limit by applying the Nystrom-type method Q 
to the Fredholm determinant E(s) = Det(I — K s ) of the 
dynamical sine kernel @, (K s -f)(x) = J Q S K(x,y)f(y)dy. 
Although this method of discretizing the Fredholm de- 
terminant by the Gaussian quadrature using m Legendre 
points {xi} and corresponding weights {wt}, 



Det(7 — K s ) ~ det [Stj — K(xi,Xj)^/wl Wjj ™ 



(3) 



FIG. 1: Ratios between Wigner-surmised LSDs P ws (s) and 
their exact forms P(s) for GOE, GUE, and GSE {fi = 1, 2, 4). 



is also a numerical approximation, uniform convergence 
to the exact LSD is proven to hold as the number m 
of quadrature points increases. The convergence is ex- 
ponentially fast, i.e., the approximation error decays as 
O(e- const m ). It is sufficient to take m to be 100 - 200, as 
the relative shifts of E(s) (for the GOE-GUE crossover) 
under the increment of m from 100 to 200 are as small 
as 10- 8 ~ 10 for £7(1), 10~ 7 ~ 8 for £7(2), 10~ 6 ~ 7 for £7(3), 
and 10~ 4 ~ 5 for £7(4). This stability ensures that numeric 
values obtained for P(s) by the Nystrom-type method at 
m = 200 (used below) can practically be regarded as ex- 
act for fitting to the actual spectra of physical systems. In 
view of this recent progress, it is now possible to examine 
the accuracies of Wigner-surmised expressions reported 
in [l| and their range of validity. In this Brief Report, I 
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FIG. 2: (Color onine) Linear and log plots of LSDs P(s) for GOE-GUE crossover evaluated by the Nystrom-type method at 
m = 200 (red liness) 0], and Wigner-surmised forms P WB (s) (blue, dotted lines) at optimally fitting values of A. A2 denotes 
the I? distance between P(s) and P WB (s). 
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FIG. 3: (Color onine) Ratios P ws (s)/P(s) for GOE-GUE 
crossover, including the two limiting cases GOE and GUE. 

shall concentrate on the simplest case of transition be- 
tween two symmetry classes, namely, LSD for the GOE- 
GUE crossover [(2.33)— (2.35) of as an example. 

Note that the crossover parameter A used in [lj (with 
si = \/tt ) is equal to \/2tt times the p parameter used 
in @, Gil- In Fig. 2, I plot LSDs P(s) at the choice of 
PI, A = 0.05, 0.1, 0.2, and 1, together with respective 
Wigner-surmised LSDs P ws (s) at A = 0.0463, 0.1828, 
0.2759, and 0.9613 that minimize the L 2 distances be- 



tween P(s) and P ws (s). These A and A 2 values are con- 
sistent with those obtained from fitting to the spectra 
of random 400 x 400 matrices (A = 0.047, 0.182, 0.276, 
0.959) shown in Fig.7 (top) of [J, which should be substi- 
tuted with Fig. 2 of this paper. The ratios P ws (s)/P(s) 
are shown in Fig. 3. 

In the range < s < 2, i.e., the validity range of the 
original Wigner surmise, the best- fitting curve P ws (s) for 
each A is hardly discernible from P(s) at corresponding 
A and serves as its good approximation, and the agree- 
ment deteriorates for s > 2. Consequently, when one 
applies Wigner-surmised LSDs to, e.g., the Dirac spectra 
of QCD-like theories at imaginary baryon chemical po- 
tential the region s > 2 should be avoided for fitting. 
And even if one restricts the fitting region to < s < 2, 
an obvious shortcoming of the Wigner surmise applied 
to the crossover ensembles lies in the difference between 
A and A: if one tries to determine the values of a or A 
in ([2]) from the spectrum of large matrices by fitting its 
LSD to the surmised expression, it inevitably involves a 
systematic error, | A/A — 1| < 0.1. In a typical example of 
extracting the pion decay constant from two-color QCD 
Dirac spectra [7|, one must instead use LSDs evaluated 
for infinite N by the Nystrom-type or analytic methods 
for fitting in order to achieve accuracy better than 10%. 
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